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I .  INTRODUCTION 


The  NEST  (N-Element  SysTem)  computer  program  is  used  to  calculate  frozen, 
equilibrium,  and  nonequilibrium  chemically  reacting,  steady,  one-dimensional 
flow.  NEST  also  can  be  used  to  calculate  the  variation  with  time  of  a  sta¬ 
tionary,  reacting  gas.  The  formulation  of  the  equations  and  numerical  methods 
used,  the  program  structure,  and  input  preparation  are  documented  in  Reference 
1.  The  Holt-Gear-Newton  implicit  method  of  integration  and  the  Emanuel-Aitken 
extrapolation  integration  method  that  are  available  in  NEST  are  described  in 
Reference  2. 

In  this  report,  the  input  preparation  and  output  options  are  documented 
and  some  sample  cases  are  presented  to  aid  the  NEST  user.  For  questions 
regarding  usage  or  chemical  species  available  on  the  thermodynamic  tape, 
contact  Karen  Foster. 


*E.  B.  Turner,  G.  Emanuel,  and  R.  L.  Wilkins,  The  NEST  Chemistry  Computer 
Program,  Vol.  I,  TR-0059(6240-20)-l,  The  Aerospace  Corporation  (30  July 
1970). 

ZW.  D.  Adams,  et  al.,  The  RESALE  Chemical  Laser  Computer  Program, 
TR-0075(5530)-5,  The  Aerospace  Corporation  (20  February  1975). 
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II.  INPUT  DESCRIPTION 


A.  DATA  STREAMS 

The  punched  card  input  for  the  NEST  computer  program  was  designed  to 
minimize  input  preparation  errors  and,  at  the  same  time,  provide  a  format  that 
is  easy  to  recall  and  prepare.  All  input  is  punched  free  form  in  columns  1 
through  70.  Columns  71  through  80  are  not  used  for  the  input  program  and  can 
contain  sequence  numbers  or  other  information. 

The  input  consists  of  many  data  streams.  A  data  stream  is  comprised  of  a 
keyword  or  keyword  followed  by  a  parameter  field.  The  keyword  identifies  a 
function  to  be  performed,  and  the  parameter  list  contains  the  details  of  the 
task,  e.g.,  the  data  stream  that  handles  all  print  options  that  can  be  con¬ 
trolled  by  input  is  identified  by  the  keyword  PRINT.  The  parameter  field  for 
this  data  stream  contains  items  such  as  the  maximum  number  of  answer  stations 
to  be  printed.  The  character  $  is  reserved  to  denote  the  start  of  a  data 
stream.  The  data  stream  can  begin  on  any  card,  in  any  column,  and  continue 
through  as  many  cards  as  necessary  for  completion  of  the  stream. 

The  program  searches  the  first  data  card  for  a  $,  which  denotes  the  be¬ 
ginning  of  the  first  data  stream.  The  search  then  continues  through  the 
card(s)  for  the  next  $,  which  denotes  the  beginning  of  the  next  data  stream. 
Then,  the  function  of  the  first  data  stream  is  determined  by  identifying  the 
keyword,  and  the  stream  is  passed  to  the  area  of  the  program  that  processes 
that  particular  function.  The  program  continues  in  this  manner  until  the  end- 
of-case  stream  is  encountered,  at  which  time  control  is  passed  to  an  area  that 
initializes  the  common  storage  area  and  calls  in  the  appropriate  computation 
overlay. 

B.  KEYWORD  RULES 

The  purpose  of  the  keyword  is  to  identify  that  data  stream.  Keywords  are 
up  to  16  alphanumeric  characters  long  and  are  chosen  to  indicate  the  function 
of  the  data  stream.  The  keyword  comprises  the  first  nonblank  characters  fol¬ 
lowing  the  $,  but  the  keyword  itself  can  contain  imbedded  blanks  to  improve 


readability.  For  example,  the  keyword  AREAVSX  can  be  punched  AREA  VS  X.  If  a 
parameter  field  follows,  the  keyword  is  terminated  by  either  a  comma  or  an 
equals  sign.  In  the  few  data  streams  that  do  not  contain  a  parameter  field, 
the  keyword  is  terminated  by  the  start  of  a  new  data  stream,  i.e.,  the  next 
$.  The  keywords  are  listed  and  described  briefly  in  Table  1. 

C.  PARAMETER-FIELD  CONSTRUCTION 

The  parameter  field  varies  for  different  types  of  data  streams.  However, 
there  are  some  general  rules  that  might  aid  in  constructing  a  parameter  field. 
Information  for  setting  up  specific  fields  is  presented  in  Section  D.  Each 
field  consists  of  a  collection  of  labels,  species  names,  and  numbers.  A  label 
is  from  I  to  16  alphanumeric  characters  long,  e.g. ,  ATM  or  ROCKET,  and  is  used 
to  specify  any  number  of  options.  Species  names  are  from  1  to  12  nonblank 
alphanumeric  characters  and  are  generally  written  in  standard  chemical  nomen¬ 
clature  (Table  2).  Thus,  gaseous  species  N2,  0,  and  F2O  are  punched  as  N2,  0, 
and  F20,  respectively.  Species  in  their  vibrational  states  such  as  ^(v  *  0) 
and  HF  (v  =*  10)  are  punched  H2(0)  and  HF(10),  respectively.  A  method  of 
punching  ionic  species  is  to  consider  a  pseudo  element  E  as  part  of  the  mole¬ 
cule  with  a  subscript  denoting  the  charge  in  terms  of  loss  or  gain  of  elec- 

+  —2 

trons.  In  this  notation,  the  three  ionic  species  O2  ,  H  ,  and  H2O  in  Table 
2  become  02*E-1,  H*E+2,  and  H20*E-3.  Condensed  species  such  as  water  and 
solid  carbon  can  be  punched  as  H20(L)  and  C  (SOLID),  respectively.  Any  non¬ 
numeric  information  within  the  parentheses  indicates  a  condensed  species,  not 
an  excited  species. 

Numbers  may  be  written  in  almost  any  way.  For  example,  the  following  are 
equivalent:  3500,  3500.,  3.5E+3,  and  +.35E4.  The  only  rule  is  that  if  a 
number  is  written  with  an  exponent,  but  with  no  decimal  point,  it  is  assumed 
that  the  decimal  point  precedes  the  first  numeric,  i.e.,  the  number  -2401E5  is 
equivalent  to  -.2401E+5. 

Blanks  in  a  parameter  field  are  ignored;  blanks  can  appear  in  species 
names,  in  numeric  fields,  and  in  labels  to  enhance  readability  without  being 
considered  part  of  the  name,  number,  or  label. 


Keyword 


Function 


Keyword 


Function 


Pg 


Pg* 


TITLE 

Problem  title 

12 

G 

Global  catalyst 

27 

FLOW 

Flow  type  description 

13 

T  RATES 

Rate  temperature 
points 

28 

ATOMS 

Atomic  composition 

15 

STEP 

Initial  step  size 

28 

QATOMS 

Pseudo  atoms 

16 

INTEG 

Integration  controls 

29 

P 

System  pressure 

16 

AREA  VS  X 

Nozzle  contour 

30 

T 

System  temperature 

17 

RADIUS  VS  X 

Nozzle  contour 

31 

RHO 

System  density 

17 

P  VS  X 

Pressure  contour 

31 

VEL 

Flow  velocity 

17 

RHO  VS  X 

Density  contour 

31 

H 

System  enthalpy 

17 

T  VS  X 

Temperature  contour 

32 

E 

System  energy 

18 

FLASH 

Flashlaap  constants 

32 

T  RANGE 

Temperature  range 

18 

I  VS  T 

Flashlamp  output 

33 

P  RANGE 

Pressure  range 

18 

PRINT 

Print  options 

34 

MOL/GM 

Composition,  mol/g 

18 

PLOT 

Plot  options 

35 

MOL  PCT 

Composition,  mol  Z 

19 

SET  GRID 

Grid  and  scaling 
controls 

39 

WT  PCT 

Composition,  wt.  Z 

19 

DUMP 

Debug  print  options 

43 

MOL/CC 

Concentration,  mol/cm^ 

19 

V  OBSERVER 

Observer  velocity 

43 

CHAMBER 

Chamber  conditions 

19 

FLAGS 

Special  controls 

43 

EXCITE 

Excited  system 

20 

END 

End  of  case 

46 

INCL 

Inclusion  of  species 

20 

NEW  CASE 

Begin  new  case 

46 

EXCL 

Exclusion  of  species 

21 

HEAT 

Heat  of  formation  update 

22 

JANAF  DATA 

Thermodynamic  data 

22 

EO 

Chemical  reaction 

24 

KF 

Forward  rates 

25 

KB 

Reverse  rates 

26 

M 

Catalyst 

26 

Table  2.  Species  Name  Input  Convention 


Species 

Method  Input 

n2 

N2 

0 

0 

f2o 

F20 

H2(v  -  0) 

H2(0) 

HF(v  =  10) 

HF(10) 

o2+ 

02*E-1 

H-2 

H2*E+2 

h2o+3 

H20*E-3 

Solid  C 

C(SOLID) 

Water 

H20  (L) 

D.  DATA  STREAM  FORMATS 

l  _ata  stream  must  begin  with  a  $  followed  by  one  of  the  keywords  listed 
in  Table  1.  If  a  parameter  field  is  required,  the  keyword  is  followed  by  a 
separator  (either  a  comma  or  an  equals  sign)  and  then  the  parameter  field. 
Items  in  the  parameter  field  are  separated  by  either  a  comma  or  an  equals 
sign. 

Data  streams  can  appear  in  any  order  in  the  input  deck,  except  that  the 
data  streams  in  the  compulsory  category  must  appear  before  all  others. 

1 .  COMPULSORY 

a.  Case  Title 

This  data  stream  contains  a  title  to  be  printed  at  the  top  of  each  page 
of  output.  The  stream  begins  with  $TITLE  ■  followed  by  a  maximum  of  72 
characters  ($  not  permitted).  Blanks  are  counted  in  the  72  character  maximum 

Example:  $ TITLE  -  HF  RUN  NO.  92 


Flow  Label 


Remarks 


SHOCK 

Nonequilibrium  shock  tube  performance  is  requested. 
This  implies  that  the  flow  area  is  constant 
(flow  area  =  1  unless  Flag  12  of  Table  14  is  input) 

EXPAND3 

Nonequilibrium  expanding  flow  specified.  A  shape 
contour,  either  area  or  radius  versus  distance, 
must  be  input.*5 

P  STREAM 

Pressure  streamtube  performance  is  specified; 
a  pressure  contour  must  be  input. 

RHO  STREAM0 

Density  streamtube  performance  is  specified; 
a  density  contour  must  be  input. 

T  STREAM 

Temperature  streamtube  is  specified;  a  temperature 
contour  must  be  input15 

Q  GAS 

The  quiescent  gas  option  applies.  This  option  must 
be  used  with  one  of  the  above. 

aNot  necessary  to  specify  throat  position. 

^The  means  of  inputting  contours  is  outlined  in  Subsection  8. 

cMust  use  Q  GAS  with  this  flow  type.  This  is  flow  type  to  use  with  flash  lamp 

Example:  $FLOW  -  RHO  STREAM,  Q  GAS 
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Table  4.  Frozen  Equilibrium  Shock  Tube  Flow  Types 


Flow  Label 

Remarks 

INCIDENT 

Calculate  frozen  and  equilibrium  conditions 
behind  the  incident  shock. 

REFLECTED 

Calculate  frozen  and  equilibrium  conditions 
behind  the  reflected  shock. 

FFI 

Calculate  the  frozen-frozen  conditions 
behind  the  incident  shock. 

FFR 

Calculate  the  frozen-frozen  conditions 
behind  the  reflected  shock. 

FF/EQ 

Calculate  frozen-frozen  incident;  frozen- 
frozen  reflected;  frozen-equilibrium 
reflected;  and  frozen-frozen,  equilibrium 
speeds,  reflected  cases  only. 

Example:  $FLOW  -  FFI,  FFR 

Frozen-frozen  incident  and  frozen-frozen  reflected  cases  are  to  be  calculated. 


Table  5.  Rocket  and  Wind  Tunnel  Flow  Types 


Flow  Label3 

Remarks 

ROCKET 

Rocket  performance  is  specified.  Use  for 
unknown  temperature.  Input  heat  of  formation. 

WIND 

Wind  tunnel  performance  is  specified. 

Use  for  pressure  and  temperature  known. 

EQ 

Calculate  equilibrium  conditions. 

FR 

Calculate  frozen  conditions. 

NONEQ 

Calculate  nonequilibrium  conditions. 

aUse  $CHAMBER  card. 


Example:  $FL0W  =■  ROCKET,  EQ,  NONEQ 

Equilibrium  and  nonequilibrium  rocket  performance  are  requested. 


Table  6.  Equilibrium  Calculations 


Flow  Label 

Remarks 

GEN  SOLVER 

General  solver,  input  a  temperature  and 
pressure  range  for  which  equilibrium 
conditions  are  to  be  found. 

MOLLIER 

Mollier  tablas  are  created  for  an  input 
temperature  and  pressure  range. 

FLAME3 

Compute  a  flame  solution  (rocket  chamber 
combustion  solution). 

VOLUME3 

Compute  constant  volume  condition. 

Input  pressure  and  density. 

aUse  $ CHAMBER  card. 

Example:  $FLOW  =  GEN  SOLVER 

c.  Atomic  Composition 

The  atomic  composition  of  the  system  must  be  entered  by  means  of  this 
data  stream.  The  keyword,  ATOMS,  is  followed  by  a  list  of  names  of  elements 
comprising  the  system.  Table  7  contains  the  complete  list  of  atoms  recognized 
by  the  program.  If  ionized  species  are  being  considered,  the  pseudo  element  E 
must  be  included  in  the  parameter  list. 

Example:  $ATOMS  =  N,  0,  AR 

The  desired  atomic  composition  consists  of  nitrogen,  oxygen,  and  argon  atoms. 


Table  7.  Element  Names  Set 


H 

NE 

K 

BR 

CS 

D 

NA 

CA 

BR/81/ 

W 

HE 

MG 

TI 

KR 

RE 

LI 

AL 

CR 

RB 

PB 

BE 

SI 

FE 

ZR 

BI 

B 

P 

CJ  (Cobalt) 

MO 

U 

C 

S 

NI 

CD 

U/235/ 

N 

'  S/34/ 

CU 

SN 

E  (electron) 

0 

CL 

ZN 

I 

F 

AR 

BR/79/ 

XE 

d.  Pseudo  Atoms 

Atoms  not  appearing  in  Table  7  can  be  included  by  means  of  the  $QATOM 
data  stream.  The  program  will  accept  as  an  element  name  the  letter  Q  followed 
by  any  letter  in  the  alphabet.  The  enthalpy,,  entropy,  and  specific  heat 
tables  for  the  new  atoms  must  be  input  by  means  of  the  $JANAF  DATA  stream. 

For  this  data  stream,  the  keyword  QATOMS  is  followed  by  the  atom  name,  the 
base  enthalpy,  and  the  atomic  weight  of  the  atom.  This  data  stream  can  be 
used  in  addition  to  the  $ATOMS  data  stream. 

Two  new  elements,  QA  and  QB,  are  defined  as  follows: 

$ QATOMS,  QA,  35,  1.5,  QB,  32.56,  15.01 
2.  STATE  VARIABLES 

a.  System  Pressure 

Keyword  P  is  followed  by  the  system  pressure  and  one  of  the  labels  ATM, 
MM,  or  PSIA.  The  cold  pressure  follows  if  $CHAMBER  is  being  used  and  mass 
flow  is  an  input  quantity. 


Example :  $P  ■  1500,  ATM 


System  pressure  is  1500  ATM. 

$P  -  837,  MM,  COLD  -  500 

System  pressure  is  837  Torr,  and  cold  pressure  is  500  Torr. 

b.  System  Temperature 

The  system  temperature,  in  degrees  Kelvin,  is  set  by  keyword  T  followed 
by  the  desired  temperature. 

Example:  $T  *  2.98E2 

c.  System  Density 

The  system  density,  in  grams  per  cubic  centimeter,  is  set  by  keyword  RHO 
followed  by  the  density.  Note  that  if  the  system  composition  is  input  in 
moles  per  cubic  centimeter,  the  program  can  be  used  to  compute  the  system 
density  and  this  data  stream  can  be  omitted. 

Example:  $RH0  *  1.9E-4 

d.  Flow  Speed 

The  flow  speed,  in  centimeters  per  second,  is  determined  by  keyword  VEL 
followed  by  the  flow  speed. 

Example:  $VEL  =  54000  • 

If  the  flow  is  one  of  the  frozen-equilibrium  shock  types  (see  Table  4),  a  set 
of  velocities  in  millimeters  per  microsecond  for  which  the  frozen-equilibrium 
conditions  are  to  be  calculated  can  follow  in  the  parameter  list. 

Example:  $VEL  ■  3,  2.5,  2,  1.5 

Conditions  are  to  be  computed  for  velocities  equal  to  3,  2.5,  2,  and  1.5. 

e.  System  Enthalpy 

The  system  enthalpy,  in  kilocalories  per  mole,  is  input  by  keyword  H 
followed  by  the  enthalpy  value. 

Example:  $H  ■  3.93 


f. 


System  Energy 

The  system  energy,  in  kilocalories  per  mole,  is  input  by  keyword  E 
followed  by  the  energy. 

Example:  $E  *  4.9 

g.  Equilibrium  Temperature  Range 

A  temperature  range  must  be  specified  for  the  general  solver  flow  type. 
This  is  set  by  keyword  T  RANGE  followed  by  an  initial  temperature,  a  final 
temperature,  and  a  temperature  decrement. 

Example: 

Answers  at  temperatures  of  3000,  2500,  2000,  and  1500  K  are  obtained  by: 
$T  RANGE  -  3000,  1500,  500 

An  answer  at  a  single  temperature  point  of  3000  requires  the  use  of 
$T  RANGE  -  3000,  3000,  0 

h.  Equilibrium  Pressure  Range 

A  pressure  range  (in  standard  atmospheres)  must  be  set  for  the  general 
solver  flow  type.  This  is  accomplished  by  keyword  P  RANGE  followed  by  an  ini¬ 
tial  pressure,  a  final  pressure,  and  a  pressure  decrement. 

Example:  $P  RANGE  ■  2.8,  1.2,  .4 

3.  COMPOSITION 

a.  Moles  per  System  Gram 

Keyword  MOL/GM  is  followed  by  sets  of  species  name  and  moles  per  gram  for 
that  species,  respectively. 

Example: 

To  input  .53  x  10”^  mol/g  of  O2,  .24  x  10”*  mol/g  of  N£,  .39  x  10”^  mol/g 

of  NO,  .252  x  10”^  mol/g  of  N,  and  .95  x  10”^  mol/g  of  0,  write 

$M0L/GM,  02  -  .53E-3,  N2  -  .24E-1,  NO  -  .39E-2,  N  -  .252E-2,  0  -  .95E-2 


b.  Mole  Percent 

Keyword  MOL  PCT  is  followed  by  sets,  each  comprising  a  species  name  and 
its  percent  amount.  The  percent  amounts  must  sum  to  100%. 

Example : 

A  composition  of  1%  water  vapor  and  99%  argon  is  specified  by 
$M0L  PCT,  H20  -  1,  AR  =■  99 

c.  Weight  Percent 

Keyword  WT  PCT  is  followed  by  sets  each  comprising  a  species  name  and  its 
weight  percent  amount.  The  percent  amounts  must  sum  to  100%. 

Example : 

A  composition  of  96.3%  fluorine  and  3.7%  carbon  monoxide  is  specified  by 
$WT  PCT,  F2  -  96.3,  CO  -  3.7 

d.  Moles  per  Cubic  Centimeter 


Keyword  MOL/CC  is  followed  by  sets,  each  comprising  a  species  name  and 
its  concentration  in  moles  per  cubic  centimeter.  Note  that  when  this  data 
stream  is  used,  the  program  can  compute  the  system  density  so  it  need  not  be 
input. 

Example :  Concentrations  of  .19  *  10~^  of  C^,  3.8  *  10-^  of  H2,  and  3.8  *  10 
of  H  are  set  by 


$M0L/CC,  CL2  -  . 19E-2,  H2  -  3.8E-6,  H  -  3.8E-9 


e.  Chamber  Conditions 


This  data  stream  is  used  to  specify  the  system  composition  for  all  rocket 
and  wind  tunnel  flow  cases,  and  the  flame  and  constant  volume  equilibrium  so¬ 
lutions.  The  keyword  CHAMBER  is  followed  by  the  chemical  name  of  the  species, 
an  optional  label  for  printout  purposes,  the  heat  of  formation  in  kilocalories 
per  mole  (not  required  for  wind  tunnel  cases),  and  the  amount  of  the  species 
being  introduced  followed  by  the  label  FUEL  if  the  species  is  a  fuel  and  the 
mixture  ratio  option  is  being  used.  In  general,  the  data  stream  is  as  follows 
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$CHAMBER,  CHEM  ■  name,  LABEL  ■  name,  HEAT  ■  number 

MASS  FLOW 
MOL/GM 

MOL  PCT  »  number,  CHEM  *  species,  etc. 

WT  PCT 

or  MIX  *  number  FUEL,  CHEM  ■  species,  etc. 

Example: 

The  system  composition  in  moles  per  gram  for  a  wind  tunnel  case  may  be 
set  by: 

$CHAMBER,  CHEM  -  02,  MOL/GM  -  .539E-3,  CHEM  -  N2,  MOL/GM  -  .241E-1, 

CHEM  -  NO,  MOL/GM  -  .394E-2 

A  rocket  case  composition  can  be  specified  in  mass  flow  (grams  per  second)  as 
follows: 

$ CHAMBER,  CHEM  -  HE,  LABEL  -  HELIUM,  MASS  FLOW  -  4.5, 

CHEM  -  SF6,  MASS  FLOW  -  2.8,  CHEM  -  02,  MASS  FLOW  -  .61 

4.  SPECIES  SPECIFICATION 

a.  Excited  Species  Inclusion 

Excited  species  are  accepted  in  the  system  if  $EXCITE  (no  parameter  field 
is  needed)  is  written. 

b.  Inclusion  of  Species 

There  are  various  ways  to  select  the  set  of  species  that  constitutes  the 
total  system  for  this  data  stream.  These  are  described  in  the  following  sub¬ 
sections. 

(1)  Case  1 

The  system  is  not  specified  by  input,  i.e.,  data  streams  $INCL  and  $EXCL 
are  not  used.  One  of  two  selection  methods  is  used  in  this  situation. 


(a)  When  an  equilibrium  or  frozen  calculation,  or  both,  is  required, 
the  program  selects  from  the  thermodynamic  data  library  tape  all  species  with 


atomic  composition  represented  by  the  set  of  elements  listed  in  the  $ATOMS  and 
$QATOMS  data  streams. 

(b)  When  only  nonequilibrium  calculations  are  required,  the  system  is 
composed  of  all  species  appearing  on  reaction  cards,  catalyst  cards,  and  com¬ 
position  cards. 

(2)  Case  2,  Inclusion  of  Species 

The  user  has  listed  the  system  species  in  a  $INCL  data  stream.  One  of 
two  selection  methods  applies. 

When  an  equilibrium  or  frozen  calculation,  or  both,  is  required,  the 
system  is  composed  of  all  species  in  the  $INCL  stream.  If  the  elements  from 
$ATOMS  are  missing  from  this  list,  the  program  will  add  them  to  the  system. 

When  only  nonequilibrium  calculations  are  required,  the  system  will  in¬ 
clude  all  species  in  the  $INCL  stream  plus  any  species  that  are  not  in  the 
$INCL  stream  but  appear  in  reaction,  catalyst,  or  composition  data  streams. 

This  data  stream  begins  with  $INCL  followed  by  species  names  separated  by 
commas  (or  equal  signs). 

Example:  $INCL  -  H,  F,  H2,  F2,  HF(0),  HF(1) 

(3)  Case  3,  Exclusion  of  Species 

Species  to  be  excluded  from  the  system  are  listed  in  the  $EXCL  data 
stream,  and  the  system  is  chosen  in  one  of  two  ways. 

When  an  equilibrium  or  frozen  calculation,  or  both,  is  required,  the 
program  begins  selecting  species  from  the  thermodynamic  data  library  as  In 
case  1(a),  but  excludes  any  species  from  the  $EXCL  data  stream.  If  one  of  the 
elements  from  $ATOMS  is  present  in  the  $EXCL  list,  its  presence  in  $EXCL  is 
ignored,  and  the  element  is  included  in  the  system. 

When  only  nonequilibrium  calculations  are  required,  the  system  Is  chosen 
as  in  case  1(b),  but  any  species  in  the  exclude  list  are  deleted.  However, 
all  species  in  reaction,  catalyst,  and  composition  data  streams  are  included 
in  the  system  even  though  they  may  appear  (assumed  incorrectly)  in  the  $EXCL 


The  methods  of  selecting  the  species  in  the  system,  as  outlined  in  the 
description  of  the  preceding  data  stream,  must  be  understood  if  this  option  is 
to  be  used.  The  data  stream  begins  with  $EXCL  followed  by  a  list  of  species 
to  be  excluded  from  the  total  system. 

Example : 

Assuming  that  0  and  C  appear  in  the  $ATOMS  list,  the  following  data 
stream  could  be  used  to  exclude  the  indicated  species: 

$EXCL  -  C02 ,  C2 ,  C3,  C5 

5.  THERMODYNAMIC  DATA 

a.  Heat  of  Formation  Update 

The  thermodynamic  library  tape  value  for  the  heat  of  formation  can  be 
changed  during  a  run  by  specifying  the  keyword  HEAT  followed  by  sets  of  a 
species  name  and  a  new  value  for  the  heat  of  formation. 

Example:  The  heat  of  formation  for  OH  and  03  can  be  set  to  9.3  and  34.2, 
respectively,  by  writing  $HEAT ,  OH  »  9.3,  03  -  34.2. 

b.  Thermodynamic  Data  Update 

This  data  stream  can  be  used  to  change  the  thermodynamic  data  from  the 
thermodynamic  data  tape  for  a  system  species  and  to  provide  data  for  a  species 
not  present  on  the  tape.  The  thermodynamics  of  any  species  comprised  of  ele¬ 
ments  defined  by  the  QATOMS  option  must  be  specified  by  this  data  stream.  The 
keywords  JANAF  DATA  are  followed  by  the  species  name,  the  heat  of  formation, 
and  a  table  containing  temperature,  molar  heat  capacity,  enthalpy,  and 
entropy.  The  temperature  values  must  be  a  continuous  subset  of  the  tempera¬ 
ture  range  on  the  tape  (see  Table  8).  If  a  range  of  temperature  values  T  is 
input  such  that  Tj  <  T  <  T2,  the  program  sets  the  thermodynamic  data  values 
for  all  T  <  Tj  equal  to  the  values  of  Tj  and  all  values  of  T  >  T£  equal  to  the 
values  of  T£« 

Example:  $ JANAF  DATA,  C2F2,  HEAT  -  5.1, 

1000,  20,  12,  78.5, 

1100,  20,  14,  80, 
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1200,  20,  15.9,  81.9, 

C2,  HEAT  =  200.0,  1000,  8.582,  6.202,  58.524 
QA02,  HEAT  -  -22.38,  1000,  23.338, 

14.412,  89.935,  1100,  23.788,  16.769,  92.181 

These  input  data  update  the  thermodynamics  for  molecules  C2F2,  C2,  and  QA02, 
as  shown  in  Table  9.  Note  that  since  only  one  temperature  was  specified  for 
C2,  all  thermodynamic  data  for  C2  are  set  equal  to  the  input  data  at  T  *  100 


Table  8.  Thermodynamic  Table  Temperature  Range 


298.15 


Table  9.  Thermodynamic  Data  Update 


Species 

C2F2 

AHF298  * 

5.1 

T 

C° 

P 

u°  -  u° 

HT  H298 

S° 

1000 

20 

12 

78.5 

1100 

20 

14 

80 

1200 

20 

15.9 

81.9 

Species 

c2 

AHF298  ’ 

200 

T 

C° 

P 

u°  _  u° 

HT  H298 

S° 

1000 

8.582 

6.202 

58.524 

Species 

qao2 

AHF298  - 

-22.38 

T 

C° 

P 

H°  -  H° 

T  298 

S° 

1000 

23.338 

14.412 

89.935 

1100 

23.788 

16.769 

92.181 

6.  REACTION  DATA 

a.  Chemical  Reaction 

The  reaction  data  stream  is  of  the  format  $EQ  XXX,  chemical  equation, 
where  EQ  XXX  is  the  keyword  portion  followed  by  a  chemical  equation.  The  term 
XXX  can  be  omitted  or  can  be  a  maximum  of  three  digits,  representing  a  unique 
serial  number  for  the  reaction.  This  serial  number  is  used  for  updating  pur¬ 
poses  only;  if  the  rate  constants  for  the  reaction  are  not  to  be  changed  in 
succeeding  cases,  the  serial  number  can  be  omitted. 

The  following  rules  must  be  observed  when  writing  a  chemical  equation. 


3* 

b. 


The  reaction  sign  (S )  is  replaced  by  an  equals  sign  (”). 

A  maximum  of  four  species  (not  including  catalysts)  appears  on  each 
side  of  the  equation. 


f. 

Examples: 


If  a  catalytic  species  is  involved,  the  catalyst  appears,  as  does  a 
species,  on  both  sides  of  the  equation,  as  MI,  where  I  is  an  inte¬ 
ger  indicating  which  catalyst  set  is  to  be  used.  In  this  case, 
catalyst  set  I  must  be  defined  by  a  catalytic  species  data  stream 
(Section  6d). 

Global  catalysts,  if  required,  also  appear  as  a  species  on  both 
sides  of  the  equation,  as  GI,  where  I  is  a  global  catalyst  set 
defined  by  a  global  data  stream  catalyst  set  (Section  6e). 

The  user  must  terminate  a  species  name  with  a  comma,  except  when 
the  species  is  the  last  species  on  either  side  of  the  reaction. 

A  species  can  appear  only  once  on  a  side. 


02  +  H2  ^  20  +  2H 

(original) 

02,  +  H2  *  20,  +  2H 

(input  version) 

M- 

02  ^  20 

(original) 

02,  +  M3  =«  20,  +  M3 

(input  version) 

In  this  example,  catalyst  set  3  must  be  defined  by  a  catalyst  data  stream. 


(3)  h2o  s:  h2o+  +  e 
H20  «  H20f,  +  E 
or  H20  -  H20*E-1 ,  +  E 


(original) 
(input  version) 


(4)  C3H603  +  Gl  *2  CH3CHO  +  CH200  +  Gl  (original) 

C3H603,  +  Gl  =  CH3CH0,  +  CH200,  +  Gl  (input  version) 

Global  catalyst  1  must  be  defined  in  a  global  catalyst  data  stream. 

(5)  H  +  H  -  2H 

This  is  unacceptable  because  it  violates  rule  f. 
b.  Forward  Rate  Constants 


»ltt] 


where  T  is  temperature  in  degrees  Kelvin  and  R  *  1.98725  (cal/mol-°K) , 

P  -  dimensionless  constant,  E  -  (cal/mol),  and  A  should  be  units  such  that 
kf  is  in  units  of  sec  and  moles/cm  .  Forward  and  backward  rate  constants  must 
be  input  and  must  immediately  follow  the  data  stream  that  defines  the  chemical 
equation.  Forward  rate  constants  are  input  by  using: 

$KF  -  Af ,  Pf,  Ef 

Examples : 

(1)  For  a  given  equation  Af  =*  .6776  x  10^,  Pf  *  -2,  and 
Ef  =-.10216  x  10^,  the  input  is: 

$KF  =  .6776E22 ,  -2,  -.10216E+6 

(2)  For  the  special  cases  in  which  temperature  is  constant  and  the  forward 
rate,  therefore,  is  constant,  only  an  Af  need  be  input. 

Thus,  if  kf  =  109,  write 

$KF  *  l.E+9 

c.  Reverse  Rate  Constants 

The  rules  that  apply  to  the  forward  rate  constants  also  apply  to  the 
reverse  rate  constants. 

Example : 

20  —1 

If  the  reverse  rate  equation  is  ■  .2  x  10  T  exp  (0/RT),  write 
$KB  -  . 2E20,  -1.,  0. 

d.  Catalytic  Species 

This  data  stream  specifies  the  species  making  up  the  catalyst  together 
with  their  weighting  factors.  The  key  word  is  MI,  where  I  is  an  integer 
indicating  this  is  the  Ith  catalytic  species  being  input.  The  parameter  field 
consists  of  sets  each  comprising  a  weighting  factor  (a  number)  followed  by  a 
species. 

Example: 

$M5  =  . 5H,  AR,  1.5F 
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defines  the  fifth  catalyst  set;  the  missing  weighting  factor  for  argon  is  as¬ 
sumed  equal  to  unity.  Note  that  $M4  =  *  defines  the  fourth  catalyst  set  to 
include  all  the  species  present  in  the  system. 

e.  Global  Catalyst  Set 

A  global  catalyst  is  defined  by 


3 

where  a  ,  br  ,  Br,  and  Cr  are  input  constants  and  Xj  =  mol/cm  of  species  j. 
The  general  fdrmat  for  this  data  stream  is 


$GR  =  B  ,  C  ,  species^  a  ,  b 

rl  ri 

species?.  a  ,  b  ,  .  ♦.,  species a  ,  b 

r2  r2  j  rj  rj 

where  r  is  an  integer  indicating  this  is  the  rth  global  catalyst.  If  Cr  *  0 
on  input,  it  is  assumed  the  entire  denominator  equals  unity  [see  example  (2) 
following] . 


Examples : 


(1)  $G1  *  1.7E14,  3.24E-3,  CL2  -  1.5,  0,  CL  =  0,  2 


implies 

G1 


<C12)1,5  (Cl)° 

1.7E14(C12)°  (Cl)2  +  3.24E-3 


(2)  $G2  *  1,  0,  CL  »  1.7,  0 

implies 

G2  -  (Cl)1,7 
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Three  temperatures,  Tp  T2,  and  T3,  to  be  used  in  the  computation  of 
either  the  forward  or  backward  rate  constant  that  has  not  been  input,  are 
specified  by: 

$T  RATES  =  Tj,  T2,  T3 

If  this  option  is  not  used,  the  program  sets  Tj  =  298  K,  T2  =  3000  K,  and 
T3  =  6000  K.  The  computation  is  performed  using  the  equation  k^/k^  =  K  where 
K  is  the  equilibrium  constant. 

Example: 

If  Tx  =  5400,  T2  =  5690,  and  T3  *  5900,  then 
$T  RATES  =  5400,  5690,  5900 
7.  INTEGRATION  CONTROLS 


a.  Integration  Step  Size  Control 

The  user  can  control  the  initial  step  size  to  be  used  in  one  of  two  ways. 

(1)  6T 

Under  this  option,  the  step  size  =  <$T/|dT/dx|.  The  quantity  in  the 
denominator  is  the  absolute  value  of  the  temperature  derivative  at  the 
starting  x.  The  program  uses  this  step  size  for  the  first  five  steps. 

(2)  Ax,  N 

Under  this  option,  both  the  initial  step  size  Ax  and  N,  the  number  of 
steps  Ax  that  is  to  be  used,  are  input. 

If  this  data  stream  is  not  used,  the  program  sets  6T  =  3  and  chooses  the 
initial  step  size  according  to  (1).  The  formats  for  this  data  stream  are: 

$STEP  =■  6T 

$STEP  =  Ax,  N 

If  two  numbers  are  input  in  the  STEP  parameter  field  it  is  assumed  that 
option  (2)  is  being  used. 


Examples : 

(1)  $STEP  *  1 
Indicates  ST  ■  1. 

(2)  $STEP  -  2.5E-8,  2 
indicates  Ax  -  2.5E-8  and  N  ■  2 

b.  Integration  Method  and  Controls 

This  data  stream  controls  the  choice  of  integration  methods.  If  this 
stream  is  not  present,  the  program  assumes  that  the  Runge  Kutta  method  is  to 
be  used  and  the  initial  step  size  is  determined  by  the  $STEP  function.  The 
data  stream  for  the  Holt-Gear  technique  is  as  follows 

SINTEG  »  HOLT,  ERROR  »  number,  HMIN  =*  number,  HMAX  *  number, 

NPRNT  ■  number,  YMIN  »  number,  EPSN  »  number 

where  ERROR  is  the  relative  accuracy  desired;  the  recommended  range  is  l.E-6  < 
ERROR  <l.E-4.  HMIN  is  the  minimum  step  size  and  HMAX  is  the  maximum  step  size 
NPRNT  is  an  optional  print  control  such  that  NPRNT  * 

0  no  extra  print 

1  print  x,  h,  y1(x),  y±' (x) 

2  print  the  above  plus  y^(x)  predicted 

3  print  the  above  plus  BND*YMAX( 1 ) ,  ICONV(I) 

4  print  the  above  during  the  Newton  iteration 

YMIN  is  an  estimate  of  absolute  zero,  y^(x)  min,  used  in  the  relative  error 
test. 

Tests  for  a  specific  problem  have  been  adequately  solved  with  the  use  of 
EPSN  ■  l.E-13;  however,  the  user  can  experiment  for  specific  problems  with  the 
use  of  smaller  or  larger  values  for  EPSN.  If  EPSN  >  l.E-10,  it  is  set  to 
l.E-13. 

Example : 


SINTEG,  HOLT,  ERROR  -  l.E-4,  HMIN  -  l.E-15, 

HMAX  -  20,  YMIN  -  l.E-11,  EPSN  -  l.E-13,  NPRNT  -  0 


The  data  stream  for  the  Emanuel-Aitken  method  is  $INTEG  =  EMANUEL,  Q1  = 
number,  HMAX  *  number,  HMIN  =  number,  DEFAULT  =  number,  RATIO  *  number,  WOBBLE 
=  number,  DOUBLE  =  number,  where  Q1  is  an  integer  such  that  if  the  integrated 
value  for  the  moles  per  gram  of  any  species  greater  than  10^*,  an  absolute  er¬ 
ror  test  is  used  for  that  species.  This  effectively  prevents  species  in  small 
concentrations  from  governing  the  step  size.  HMAX  and  HMIN  are  the  maximum 
and  minimum  allowable  step  sizes,  respectively.  The  error  is  kept  between  a 
value  equal  to  DEFAULT  and  DEFAULT/RATIO.  The  variance  in  the  estimate  of  the 
derivatives  is  kept  less  than  a  value  equal  to  WOBBLE,  and  the  doubling 
criterion  must  be  passed  DOUBLE  steps  before  the  step  size  is  increased. 

Keywords  DEFAULT,  RATIO,  WOBBLE,  and  DOUBLE  need  not  be  specified.  The 
program  default  values  are  shown  on  the  example  card. 

Example : 

$INTEG  -  EMANUEL,  Q1  =■  -8,  HMAX  -  2,  HMIN  =  l.E-4,  DEFAULT  =  l.E-4, 

RATIO  =  20,  WOBBLE  =*  .03,  DOUBLE  =*  2 

8.  CONTOURS 

For  each  of  the  five  contours  described  in  the  following  paragraphs,  a 
set  of  coordinates  describing  the  contours  is  input.  A  maximum  of  100  coor¬ 
dinates  can  be  input.  If  only  one  point  is  input,  the  contour  is  assumed 
constant.  If  two  points  are  input,  the  program  interpolates  linearly  to  find 
intermediate  results.  Quadratic  interpolation  is  used  when  more  than  two 
points  are  input. 

a.  Area  Contour 

When  this  data  stream  is  used,  it  is  assumed  that  the  flow  is  a  rocket  or 
wind  tunnel  nozzle  or  an  expanding  flow  type.  The  area  contour  (in  square 
centimeters)  is  input  as: 

$AREA  VS  X  *  Aj,  x^,  A2,  X2«*»*» 

where  A^  is  the  contour  area  at  distance  x^.  (x^  is  the  time  if  the  quiescent 

gas  option  is  on.) 


Example : 

$AREA  VS  X  *  1,  0,  5.4,  1,  45,  10  specifies  area  values  1,  5.4,  and  45  at 
distances  0,  1,  and  10. 

b.  Radius  Contour 

This  option  can  be  used  instead  of  the  area  contour  whenever  it  is  more 
convenient  to  input  the  radius  contour.  The  radius  contour  (in  centimeters) 
is  input  as: 

$ RADIUS  VS  X  ■  ,  x^  ,  r2  »  *2*  *  *  * 

where  r^  is  the  value  of  the  radius  at  distance  x^. 

c.  Pressure  Contour 

When  this  data  stream  is  used,  it  is  assumed  that  the  flow  type  is  a 
pressure  streamtube  (P  STREAM).  The  pressure  contour  (in  standard  atmo¬ 
spheres)  is  input  as: 

$P  VS  X  m  pp  Xj,  P2 »  x2 , . « • 

where  p^  is  the  system  pressure  at  distance  x^.  (x^  is  time  if  the  quiescent 

gas  option  is  on.) 

Example : 

$P  VS  X  ■  .38,  0,  .32,  4.E-4,  .11,  7.E-3  specifies  pressures  of  0.38, 
0.32,  and  0.11  atm  at  x  ■  0,  4.E-4,  and  7.E-3,  respectively. 

d.  Density  Contour 

When  this  data  stream  is  used,  it  is  assumed  that  the  flow  type  is  a 
density  streamtube  (RHO  STREAM).  The  density  contour  (grams  per  cubic 
centimeter)  is  input  as: 

$RH0  VS  X  ■  p  ^ ,  P2 ’  X2 *  #  *  * 
where  is  the  system  density  at  distance  (time)  x^ 

Example : 

$RH0  VS  X  -  1.09E-5,  0,  1.08E-5,  10 

specifies  densities  1.09E-5  and  1.08E-5  at  points  0  and  10,  respectively. 
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e.  Temperature  Contour 

When  this  data  stream  is  used,  it  is  assumed  that  the  flow  type  is  a 
temperature  streamtube  (T  STREAM).  The  temperature  contour  (in  degrees 
Kelvin)  is  input  as: 


$T  VS  X  -  Tj_,  xj^,  T2,  x2  ••• 
where  T^  is  the  temperature  at  distance  (time)  x^. 
Example: 


$T  VS  X  ■  3000,  0  indicates  the  temperature  is  3000  K  at  x  ■  0,  and  since 
only  one  temperature  point  is  input,  the  temperature  is  assumed  to  be  constant. 

9.  FLASHLAMP  OPTIONS 

a.  Photodissociation  Reactions 

Equation  IIB-15  of  Ref.  1  was  modified  as  follows: 


ph 


ClZi[I(t) ] 
P 


1 


exp 


> 


This  change  was  made  to  permit  more  accurate  modeling  of  conditions  wherein  a 
species  in  an  excited  state  is  formed,  rather  than  dissociation  products.  De¬ 
fault  values  of  ga  ■  gb  ■  1  and  nb  -  0  are  built  into  the  program. 


A  maximum  of  five  photodissociation  reactions  can  be  input.  The  format 
of  the  photodissociation  reactions  data  stream  is  $FLASH,  ABSORB,  Z1  ■  number, 
Z2  *  number,  VA  ■  number,  GA  ■  number,  GB  «  number,  NB  ■  species  name,  species 
name  -  c^,  species  name  “  c2,  ABSORB,  etc.  The  order  following  ABSORB  may 
vary. 


ABSORB 

Z1 

Z2 

VA 


label  Indicating  the  start  of  a  parameter  string  containing 
constants  for  one  reaction 

maximum  number  of  moles  of  laser  gas  available  in  moles  per 
second  per  cubic  centimeter 

laser  gas  absorption  parameter  in  cubic  centimeters  per  mole 
photon  frequency  in  sec”^ 
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GA  -  degeneracy  of  the  state,  ci  -  -1  (dimensionless). 

GB  ■  degeneracy  of  the  state,  c^  ■  +1  (dimensionless). 

Default  value  ■  1. 

NB  ■  name  of  species  produced  by  absorption.  Default  value  ■  0. 

The  sets  of  pairs  of  species  name^  *  c^  contain  the  name  of  a  species  and 
its  value  of  c^  where  c^  is: 

-1  if  the  species  is  dissociated  by  photons 

4-1  if  the  species  is  produced  by  photons 

+2  if  2  molecules  of  species  1  are  produced  by  dissociating  one  molecule 
(For  all  other  species,  no  is  input.) 

Example: 

This  example  Illustrates  the  input  for  a  case  in  which  only  one  reaction 
is  needed;  F2O  is  the  absorbing  species  and  no  species  b  is  needed. 

$ FLASH,  ABSORB,  Z1  -  .129,  Z2  -  2.3E5, 

VA  -  1.55E15,  F  ■  1,  FO  *  1,  F20  -  -1 

This  example  illustrates  a  case  in  which  only  one  reaction  is  needed, 
HBR/79/(4)  is  the  absorbing  species. 

$ FLASH,  ABSORB,  Z1  -  .15,  Z2  -  2.29E6,  VA  -  7.4E8, 

HBR/79/ (4)  -  -1,  HBR/79/( 13)  -  +1,  NB  -  HBR/79/(13), 

GA  -  1,  GB  -  .778 

b.  Flash lamp  Output  Function 

A  table  of  the  normalized  dimensionless  flashlamp  output  (I)  as  a 
function  of  time  (t)  must  be  input  for  each  photodissociation  reaction  input 
in  the  FLASH  data  stream.  A  maximum  of  100  points  can  be  input  in  each  table 
as  follows: 

$1  VS  T  ■  1^,  t^ ,  I2 ,  t2 , « • « 

For  example,  if  two  photodissociation  reactions  are  input,  then  two  separate 
output  tables  must  be  input  in  the  same  order  that  the  reactions  appear  on  the 
$FLASH  data  stream. 


-  J 

•'  ) 

I 


* 


- 1 

.-■1 


. 


ft* 


L 
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ft; 


$1  VS  T  -  0,  0,  .2,  2.5E-6,  .4,  .5E-6,  .6,  7.5E-6,  .78,  l.E-5,  .9,  1.2E-5 


$1  VS  T  -  0,  0,  .3,  -  3.0E-6,  .5  6.2E-6,  .6,  7.6E-6,  .8,  1.1E-5,  .1, 
1.3E-5 

The  information  in  the  following  two  sets  of  tables  is  implied. 


T(  1)  t(D 

0  0 

.2  2.5E-6 

.4  5.0E-6 

.6  7.5E-6 

.78  1.0E-5 

.9  1.2E-5 

10.  OUTPUT  OPTIONS 

a.  Printed  Output  Gontrols 

This  data  stream  contains  print  control  options  and  is  of  the  form 

SPRINT  -  THERM,  RATES,  SUMMARY,  INTERVAL  -  number,  MAX  ANSWERS  -  number, 

X  STOP  ■  number,  THROAT 

The  order  in  which  these  options  appear  can  vary. 

THERM 

If  this  label  is  present,  it  causes  a  printout  of  the  thermodynamic  data 
library  tables  for  each  species  in  the  system. 

RATES 

If  this  label  is  present,  the  values  of  KF,  KB,  LF,  LB,  and  (LF  -  LB)  are 
printed  at  each  answer  station. 


jC  2)  t(2) 

0  0 

.3  3.0E-6 

.5  6.2E-6 

.6  7.6E-6 

.8  1.1E-5 

1.  1.3E-5 


SUMMARY 


If  Chls  label  Is  present,  the  engineering  summary  tables  are  printed  for 
rocket  and  wind  tunnel  performance  only. 

INTERVAL  ■  number 

This  label  specifies  the  Integration  step  where  answers  are  to  be 
printed. 

INTERVAL  *  I  causes  printout  at  the  end  of  each  integration  step. 

MAX  ANSWERS  -  number 

This  label  specifies  the  maximum  number  of  answer  sets  to  be  printed. 

X  STOP  *  number 

This  indicates  the  last  answer  station,  i.e.,  value  of  the  independent 
variable  (distance  or  time). 

THROAT 

This  label  can  be  used  for  rocket  and  wind  tunnel  runs  for  answer 
printout  at  the  chamber  and  throat  locations  only. 

Example : 

$ PRINT,  MAX  ANSWERS  -  100,  INTERVAL  -  20,  RATES,  X  STOP  -  45 
b.  Plot  Output  Control 

Plots  can  be  obtained  for  the  flow  types  listed  in  Table  3.  The  keyword 
PLOT  is  followed  by  the  list  of  variables  to  be  plotted.  The  individual  vari¬ 
ables  must  be  separated  by  either  a  comma  or  a  slash,  a  comma  indicates  that 
the  next  variable  is  to  be  plotted  on  the  same  set  of  axes  as  the  preceding 
variable  (a  maximum  of  nine  variables  can  be  plotted  on  a  single  set  of  axes), 
and  a  slash  indicates  that  the  following  variable  is  to  be  plotted  on  a  new 
set  of  axes. 

The  designations  used  to  identify  the  variables  are  derived  from  the 
printed  output.  The  basic  designations  recognized  by  the  program  are  shown  in 
Table  10. 


p 

Pressure 

T 

Temperature 

H 

Enthalpy 

S 

Entropy 

R 

Density 

W 

Weight  flow 

CP 

Specific  heat  at  constant 
pressure 

G 

Ratio  of  specific  heats 

M 

Mach  number 

V 

Flow  velocity 

AREA 

SIGMA1 

GAM1 1 

GAM21 

LAMBDA l1 

LAMBDA2 

Flow  area 

A 

Acoustic  velocity 

I 

Normalized  flashlamp 
output  function 

The  molar  concentration  of  a  species  is  designated  by  the  name  of  the 
species,  which  must  be  identical  to  the  name  as  it  appears  on  the  thermo¬ 
dynamic  data  tape.  An  alternative  designation  is  "CONC(Sp  S£ »•..)"  where  the 
Sp  S2»  etc.,  are  species  names.  The  latter  form  must  be  used  for  the  species 
H,  S,  and  I  to  avoid  confusion  between  these  species  and  the  designations  for 
enthalpy,  entropy,  and  the  normalized  flashlamp  output  function. 

Partial  derivatives  can  also  be  plotted;  their  proper  designations  are 
listed  in  Table  11. 

Note  that  the  partials  for  temperature,  area,  density,  and  pressure  can 
be  combined,  e.g.,  D(R,T)/DX  will  cause  the  partial  derivatives  of  both 
density  and  temperature  to  be  plotted. 

To  plot  the  rate  of  production  of  a  chemical  species,  the  name  of  the 
species  is  inserted  within  the  parentheses  of  the  designation  D(  )/DX.  For 
flashlamp  cases,  D(  )/DX  or  CHEM  D(  )/DX  designates  photon  rates  of  produc¬ 
tion.  More  than  one  species  name  can  be  included  within  the  parentheses. 

To  plot  rate  equation  data,  the  appropriate  designation  is  used  with  the 
desired  equation  numbers  inserted  within  the  parentheses.  The  designations 
recognized  by  the  program  are  listed  in  Table  12. 

Example; 

$PLOT,  KF(3, 10/1) ,  LF ( 5 ) / H ,  S/H2,  CONC(H)/D(HF(0) ,  F0)/DX/AREA 
will  produce  six  separate  plots,  as  follows: 

Plot  number  Contents 

1  KF  rate  data  for  equations  3  and  10 

2  KF  rate  data  for  equation  1  and  LF  rate 

data  for  equation  5 

3  Enthalpy  and  entropy 

4  Molar  concentrations  for  H2  and  H 

5  Rates  of  production  for  HF(0)  and  F0 

6  Flow  area 


Table  11.  Derivative  Designations 


U 


Plot  Designation 

Variable 

■ 

D(T)/DX 

Partial  derivative 

of 

temperature  J 

D(A)/DX 

Partial  derivative 

of 

area 

D(R)/DX 

Partial  derivative 

of 

density 

D(P)/DX 

Partial  derivative 

of 

pressure 

Lj 

DX 

Step  size 

« 

1 

• 

i-  . 
r.  ■  . 

DTIME 

Step  size  for  quiescent  gas 

Table  12.  Rate  Equation  Designations 


a.  The  maximum  number  of  items  that  can  be  plotted  is  100 

b.  If  $SETGRID  is  not  used,  the  plot  program  automatically  selects  the 
grid  and  scaling.  If  more  than  one  item  is  to  be  plotted  on  the 
same  set  of  axes,  the  grid  and  scaling  will  be  selected  on  the  basis 
of  the  first  item  only. 

c.  All  data  points  used  in  the  plots  are  automatically  taken  from  the 
normal  program  output. 

d.  Excited  species  are  permitted.  The  species  names  are  employed  in 
the  same  manner  as  unexcited  species  names. 

e.  The  program  will  not  interpret  the  slash  mark  that  appears  as  a  part 
of  the  partial  derivative  and  rate  of  production  designations  as 
signifying  that  a  new  plot  will  be  started. 

f.  Mo  comma  or  slash  should  be  placed  after  the  last  item. 

g.  All  variables  are  plotted  as  a  function  of  particle  time. 


c.  Plot  Grid  and  Scaling  Controls 


This  data  stream  can  be  used  to  override  the  automatic  scaling  of  axis 
and  the  choice  of  grid  by  the  NEST  plot  program.  The  keyword  SETGRID  is 
followed  by  the  frame  number,  the  grid  type,  and  scaling  information  as 
follows: 


$SETGRID,  FRAME  -  1,  GRID  -  XXX,  XMIN  -  value,  XMAX  -  value,  XDELTA 


value,  YMIN  -  value,  YDELTA  ■  value,  YSCALE  ■  value,  FRAME 


GRID 


The  GRID  patterns  that  can  be  used  are  listed  in  Table  13. 


XMIN 

XMAX 

XDELTA 


starting  point  for  X  axis  of  plot, 
will  be  discarded. 


Any  value  less  than  XMIN 


maximum  value  along  the  X  axis.  Any  value  greater  than  this 
will  be  discarded. 


number  of  units  per  inch  along  the  X  axis.  If  optimum  scaling 
along  the  X  axis  is  desired,  then  appropriate  values  for  XDELTA 
are  -1  for  10-in.  plots,  -2  for  20-in.  plots,  -3  for  30-  ,. 

plots,  and  -4  for  80-in.  plots 


YMIN 

YDELTA 

YSCALE 


lowest  or  starting  point  for  the  Y  axis, 
number  of  units  per  inch  along  the  Y  axis. 


for  log  plots  only,  minus  the  number  of  inches  per  cycle.  See 
Table  13  and  use  YDELTA  -  0. 


Table  13.  Plot  Grid  Selection 


Type 


Grid  Number 


Description 


No  grid 

300 

Square  grids 

301 

X  10  divisions  per  inch 

Y  10  divisions  per  inch 

302 

X  20  divisions  per  inch 

Y  20  divisions  per  inch 

3A1 

1-in.  border  on  left 

X  10  divisions  per  inch 

Y  10  divisions  per  inch 

3A2 

1-in.  border  on  left 

X  20  divisions  per  inch 

Y  20  divisions  per  inch 

3B1 

1-in.  border  on  left  and  right 
X  10  divisions  per  inch 

Y  10  divisions  per  inch 

3B2 

1-in.  border  on  left  and  right 
X  20  divisions  per  inch 

Y  20  divisions  per  inch 

3CM 

1-in.  border  on  left,  0.25-in. 

border  on  right 

(grid  size  17  x  9.84  in.) 

X  10  divisions  per  0.79  in. 

(5  divisions  per  centimeter) 
Y  10  divisions  per  0.79  in. 

(5  divisions  per  centimeter) 

Rectangular  grids 

3AS 

1-in.  border  on  left,  0.7-in. 
border  on  right 

X  5  divisions  per  0.85  in. 

Y  10  divisions  per  inch 

3BS 

1-in.  border  on  left 

X  10  divisions  per  0.8  in. 

Y  10  divisions  per  inch 

Semi-log  grids 

3AL 

1-in.  border  on  left  and  right 
X  10  divisions  per  0.6  in. 

Y  4  cycles,  2.5  in.  per  cycle 

40 


Table  13.  Plot  Grid  Selection  (con't.) 


3BL 

1-in.  border  on  left  (grid 
size  17  x  9.8  in.) 

X  10  divisions  per  inch 

Y  7  cycles,  1.4  in.  per  cycle 

3DL 

1-in.  border  on  left  and  right 
X  10  divisions  per  inch 

Y  5  cycles,  2  in.  per  cycle 

3KL 

1-in.  border  on  left 

X  20  divisions  per  inch 

Y  3  cycles,  3.33  in.  per  cycle 

3ML 

1-in.  border  on  left 

X  20  divisions  per  inch 

Y  8  cycles,  1.25  in.  per  cycle 

3NL 

1-in.  border  on  left  and  right 
X  2  cycles,  7.5  in.  per  cycle 

Y  10  divisions  per  inch 

30L 

1-in.  border  on  left  and  right 
X  3  cycles,  5  in.  per  cycle 

Y  10  divisions  per  inch 

3PL 

1-in.  border  on  left 

X  20  divisions  per  inch 

Y  5  cycles,  2  in.  per  cycle 

Log-log  grids 

3CL 

1-in.  border  on  left  and  right 
X  2  cycles,  7.5  in.  per  cycle 

Y  2  cycles,  5  in.  per  cycle 

3EL 

1-in.  border  on  left 

X  8  cycles,  2  in.  per  cycle 

Y  5  cycles,  2  in.  per  cycle 

3FL 

2-in.  border  on  left,  1-in. 
border  on  right 

X  7  cycles,  2  in.  per  cycle 

Y  5  cycles,  2  in.  per  cycle 

3GL 

1-in.  border  on  left  and  right 
X  5  cycles,  3  in.  per  cycle 

Y  5  cycles,  2  in.  per  cycle 

Table  13.  Plot  Grid  Selection  (con't.) 

3HL  2-in.  border  on  left 

X  5  cycles,  3  in.  per  cycle 

Y  5  cycles,  2  in.  per  cycle 

3JL  1-in.  border  on  left  and  right 

X  2  cycles,  7.5  in.  per  cycle 

Y  5  cycles,  2  in.  per  cycle 


Example : 


$SETGRID,  FRAME  -  1,  GRID  -  3BL,  XDELTA  -  -1, 

YMIN  =  l.E-14,  YDELTA  =  0,  YSCALE  =  -1.4,  FRAME  -  2 

implies  two  plots  are  to  be  made  with  automatic  X  axis  scaling  on  7-cycle 
semi-log  paper  with  a  minimum  of  1  x  10“^  for  the  Y  axis. 

Dimensions  for  all  grids,  including  margins,  are  17  in.  (x)  by  10  in. 

(y),  except  where  otherwise  noted,  with  5/8  reduction  for  CDC  HARDCPY. 

d.  Program  Trace 

The  program  trace  flag,  ID3R,  can  be  set  by  the  following: 

a.  $D(JMP  =■  ON 
sets  ID3R  -  -1 

b.  $DUMP  -  OFF 
sets  ID3R  *  0 

c.  $DUMP  *  ON,  N 
sets  ID3R  »  N 

e.  Observer  Velocity 

The  velocity  relative  to  the  observer  (in  centimeters  per  second)  is  used 
to  compute  and  print  particle  time  for  Table  3  flow  types.  The  input  format 
is 

$V  OBSERVER  -  number 
Example:  $V  OBSERVER  -  25000 
11.  SPECIAL  PROGRAM  CONTROLS 
a.  Program  Flags 

This  option  permits  the  user  to  alter  some  of  the  built-in  flags  in  the 
program.  Table  14  contains  a  list  of  these  flags,  which  are  input  as  follows: 

$FLAGS ,  f  ^  *  0^  , f 2  9  n2  ,  •  •  • 

f * ,  f-»  are  flag  numbers  that  are  to  be  set  to  n, ,  nj,  etc. 


Table  14.  Special  Option  Flags 


Flag 


Description 


2  Starting  temperature  for  flame  solution  if  different  from  3000  K. 

2 

11  Rocket  exit  area  (cm  ). 

12  Value  of  inlet  area  for  nonequilibriura  runs  (not  including  rocket 
and  wind  tunnel);  default  value  is  1  cm  . 

13  Internal  energy  computed  from  input  conditions  is  increased  by  an 
amount  equal  to  the  value  of  flag  13  before  the  constant  volume 
combustion  calculation  is  performed. 

14  If  flag  14  is  nonzero,  the  L*r  and  arrays  are  computed  based  on 
the  equilibrium  chamber  conditions  in  the  rocket  overlay.  These 
values  are  printed  and  the  program  immediately  returns  to  the  input 
overlay. 

15  Stability  parameter  B,  set  ■  1  unless  flag  15  is  used. 

17  If  flag  17  is  nonzero,  extra  printout  is  provided  during  the 
calculation,  in  the  input  overlay,  of  the  state  variables  not 
specified  on  input  cards. 

18  Value  of  the  mass  flow  in  rocket  and  wind  tunnel  cases.  Alternative 
use:  if  *0,  return  to  the  input  overlay  for  the  next  case  after 
detecting  errors  in  ESHOCK  (otherwise  exit). 

2 

19  Contains  the  throat  area  (cm  )  for  rocket  and  wind  tunnel  cases. 
Alternative  use:  if  *0,  punch  concentrations  of  excited  species  in 
laser  constant  density  runs. 

20  *0,  print  L^r  and  arrays  at  every  step  during  the  ROCKET  overlay 
equilibrium  computations. 

29  Normally  *  100;  if  the  max/min  ratio  of  a  variable  to  be  plotted  is 
greater  than  the  input  value  of  flag  29,  use  log  plots. 

30  If  the  system  freezes  in  ESHOCK,  i.e.,  if  the  percent  temperature 
difference  between  corresponding  frozen  and  equilibrium  solutions  is 
less  than  flag  30,  then  compute  only  the  frozen  conditions  for  the 
remaining  velocities. 

31  The  computation  is  terminated  in  SHOCK  when  the  Mach  number  is  less 
than  flag  31. 
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32  In  the  rocket  and  wind  tunnel  summary  table,  flag  32  (normal  value  = 

1)  is  the  amount  the  area  ratio  is  incremented  to  print  the  next  set 
of  answers. 

34  Any  area  downstream  of  the  throat  (set  equal  to  3  normally).  See 
flag  35. 

35  Number  of  print  stations  to  be  used  between  the  throat  and  the  area 
of  flag  34  (set  equal  to  32  normally). 

37  *0,  for  a  concise  trace  printout  in  ROCKET  and  NEQUIL. 

38  If  nonzero,  flag  38  equals  the  x  coordinate  at  which  the  card  type 
10  answer  interval  is  to  be  overridden  and  printing  begun  of  answers 
at  every  step  in  overlays  SHOCK  and  CALC. 

39  In  SHOCK  and  CALC,  flag  39  equals  the  x  coordinate  at  which  the  dump 
flag  ID3R  is  to  be  set  equal  to  -1.  In  ESHOCK,  flag  39  equals  the 
decimal  number  R.C,  where  R  represents  the  run  number  and  C  the  case 
number  at  which  the  dump  flag  ID3R  is  set  equal  to  -1. 

40  If  6  (divergence  point  parameter)  is  less  than  flag  40,  departure 
from  equilibrium  has  occurred  in  ROCKET.  (Flag  40  is  normally  set 
equal  to  0.015). 

42  Flag  42  equals  the  maximum  allowable  integration  step  size  (normally 

set  equal  to  0.2  cm). 


9 


Maximum  number  of  NEST  printout  pages  permitted 


Example : 


$FLAGS,  42-5 

sets  the  maximum  allowable  integration  step  size,  i.e.,  flag  42  equal  to  0.5  cm. 
12.  EDIT  CONTROLS 

a.  End  of  Case 

Each  case  must  end  with  $END. 

b.  New  Case 

The  new  case  data  stream  is 
$NEW  CASE 

All  complete  new  cases,  with  the  exception  of  the  first  case,  must 
begin  with  this  data  stream. 

c.  Multiple  Case  Runs 

Multiple  case  runs  can  be  handled  by  the  input  program.  The  run  may 
consist  of  the  stacking  of  complete  cases,  the  updating  of  previously  input 
cases,  or  both;  however,  the  deck  setup  and  assumptions  differ  for  the  two 
situations. 

(1)  Complete  Case  Stacking 

A  complete  case  consists  of  all  the  necessary  input  data  for  the  solution 
of  one  set  of  conditions  for  one  flow  type.  The  method  of  stacking  complete 
cases  is  illustrated  in  Fig.  1.  Note  that  all  cases  but  the  first  begin  with 


COMPLETE  SET  OF  DATA 
FOR  CASE  1,  ENDING 
WITH  $END 


$NEW  CASE  FOLLOWED 
BY  COMPLETE  SET  OF 
DATA  FOR  CASE  2, 
ENDING  WITH  $END 


$NEW  CASE  FOLLOWED 
BY  COMPLETE  SET  OF 
DATA  FOR  LAST  CASE, 
ENDING  WITH  $END 


Complete  Case  Stacking 


In  Chis  situation,  the  complete  set  of  data  for  the  first  case 
(terminated  by  $END)  is  followed  by  sets  of  items  (each  set  terminated  by 
$END)  to  be  updated  in  the  previous  case.  The  quantities  that  can  be  updated 
are  restricted;  the  following  data  streams  cannot  be  repeated  after  the  first 
case  is  input. 

a.  $FLOW 

b.  $ATOMS 

c.  $QATOMS 

d.  $EXCITE 

e.  $INCL 

f.  $EXCL 

g.  $HEAT 

h.  $JANAF  DATA 

i.  $EQ 

j.  $M 

k.  $G 

The  major  items  in  the  system,  i.e.,  atomic  composition,  flow  type, 
system  species,  and  reaction  equations,  therefore,  cannot  be  updated.  In  Fig 
2  the  method  of  using  the  update  feature  is  illustrated.  The  first  complete 
case  is  input,  terminated  by  $END,  and  followed  by  a  set  of  data  updating  the 
first  case. 

In  the  second  case  to  be  executed,  the  title,  system  temperature,  and 
forward  rate  for  the  reaction  with  serial  number  15  are  updated. 


$TITLE  *  CASE  1 


$  ATOMS  *  H,  F,  AR 

$FLOW  *  RHO  STREAM,  $T  *  5000 

$EQ  15,  2H  -  H2  $KF  3  5E19,  -1,  0 


$END 


III.  CONTROL  CARDS 


The  following  control  cards  must  be  used  to  run  the  NEST  Program 

ATTACH,  TAPE5 ,  8TAPE5 ,  ID  -  09680. 

ATTACH,  MAIN,  4NESTMAIN,  ID  -  06629. 

ATTACH,  L,  4NESTSUBS,  ID  »  06629. 

ATTACH,  PLOTLB,  3FTNPL0TLB. 

LDSET,  LIB  -  L/PLOTLB,  PRESET  »  ZERO,  MAP  *  0. 


IV.  SAMPLE  INPUT  SHEETS 


Seven  sample  cases  are  presented  in  this  section.  These  cases  were 
prepared  for  demonstration  purposes  only;  it  should  not  be  assumed  that  any 
quantities  such  as  rate  constants  or  heats  of  formation  are  accurate. 


The  first  sample  input  sheet  is  a  typical  nonequilibrium  problem.  The 
flow  type,  specified  on  the  second  card,  Is  a  constant  pressure  case.  Since 
the  composition  is  input  in  mol/cm  and  the  temperature  is  specified  on  card 
3,  the  program  can  compute  the  pressure  of  the  system.  This  pressure  infor¬ 
mation  completes  the  program  requirement  that  two  state  variables  be  input, 
and  since  the  pressure  is  not  specified  in  the  input,  it  is  assumed  that  this 
is  a  constant  pressure  problem.  Ihe  $INTEG  data  stream  was  omitted,  so  the 
default  Runge-Kutta  integration  method  will  be  used  with  the  initial  6T 
specified  on  the  $STEP  data  stream.  Card  5  sets  the  output  controls  for 
printout  at  every  second  integration  step;  the  problem  is  completed  when 
either  the  variable  X  exceeds  10  cm  (i.e.,  XSTOP)  or  a  maximum  of  100  answers 
is  printed. 


CARO 

NO. 

iT:TLt=F-vE3SURR  S  r  RE  A  M  TU^t 

CAKO 

NO. 

2 

i  AT0MS=C,0,H,N  .r-D*  =  P  STREAM 

CARO 

NO  . 

3 

wT-513  7£L=2.9?531£S 

CARD 

NO. 

t 

i  -S  t-p= .  i.  *  j  i 

CAF.O 

NO. 

5 

RINT.INTE  FVA_=2,  UX  ANSWE-<S=1JC,XST  JP  =  i: 

CARO 

NO. 

6 

;N0t/CC,C02=.53t7£-9t H£= . i 9 9 3£- 9, H= . l lu£- 1 C 

> 

CARO 

NO. 

7 

H2d  =  .93eE-9-,0H=.67E-ll,G  =  .  2s7I-12,C0=.6CaE- 

9 

CARO 

NO. 

3 

vHlOT.T/COMC  (H.O.OH,  H20.  HO  2  ) /CO  .  CO  ?,  H  3,  M  2 . 0 
jSETGrI0»P^AMEs1,"?I3  =  3GZ.X 9  J  N= 3 , XM AX  =  •  C  2» 
XQElT  A  =  . 3.2,YHIN=5JJ, YGELT 4=5 . FRA N£=2 , G R 19= 

** 

c 

CARO 

NO. 

9 

CARD 

NO. 

i; 

331 , 

CAC  0 

NO. 

li 

Y'-!lN=i.E-lt,Y0£LTA  =  3,  YSCAL£  =  -i.t , FRAMES 3 

CARO 

NO. 

12 

iHl=H.N,0,HO|H20,H2,00,C02fNC,N2,02 

CARO 

NG. 

13 

ic0=C02, *H1=C0, *5,  *11  iK6  =  . l£l?, 3 25 t 3 

• 

CARO 

NO. 

It 

c  *0=H20  ,*.'1t=GH,  +3, *11  ;<d=.  IE 20,-1. 

CARD 

NO. 

15 

it 1=H2,*M1= 2H,*Ml  tKd=. 75 £19, -1 . 

CARO 

NO. 

16 

;,£Q=C02,  *H  =  C0  ,*0H  ;<9=3.1£li,  0,  -6  CC 

• 

CACQ 

NO. 

17 

EQ=H 20,  *H=QH,*H2  *K3=.  6E12,  .  5, -5  *3 

L  • 

CARO 

NO. 

13 

i£Q  =  H0  2, *H=M2,* 0  2  ; K F= 1 . Ell , . 5,  -6  £3 

r 

*  • 

CARO 

NO. 

19 

;EQ=02,*H2=2 OH  r  = . IF  15 , C ,- . 7 £5 

CARO 

NO. 

23 

f£0=N2,*02=2N0  *<F=. 2E13,:,-. 122 

E  6 

CARO 

NO. 

21 

ic.3=C0,  *02  =  0,*: 02  jKF=.  33E13,  J,  -  .  51 

£5 

CARO 

NO. 

22 

i  END 

uno 


All  eight  cases  described  in  Table  4  are  to  be  calculated.  Either  the 
system  pressure  or  density  must  be  input  and  the  ambient  temperature  must  be 
specified;  in  this  case  T  =  299.2  K.  The  eight  cases  are  to  be  calculated  for 


each  input  velocity. 


t- 

. 


ARD  NO. 
ARD  NO. 
ARO  NO. 
CARO  NO. 
CARO  NO. 


1 

7 

3 

4 
3 


:TITLr_-=r:QUlLlFKl'JM/-  =>CZEN  SH 
;.FlJW  =  InC1  It  NT,  -  E-  utCTEO  i 

:R*.G65739-»,ATM  ST=299.2 

iisCL  i=CT  .H20  =  l.  A-  =99 


SHOCK  TU IE 
‘ AT  C  M  S=H , 0 t A 


*5 

•  c. 


This  is  an  equilibrium  rocket  performance  problem  and  is  so  specified  on 
card  2.  The  system  pressure  must  be  input  for  a  rocket  case.  The  heats  of 
formation,  used  to  calculate  the  system  enthalpy,  are  input  for  the  components 
that  appear  on  the  CHAMBER  card.  Card  5  contains  a  very  complex  chemical  name 

for  one  of  the  system  species,  c. i 17^9 . 5855°1 . 684N. 1384*  This  is  the  only 
data  stream  where  fractional  subscripts  may  appear  in  the  species  name.  Card 
3  sets  flag  11,  the  nozzle  exit  area,  equal  to  50. 


Cft-0  NO.  i  $  T;TLr:=EQUlLl"RlJ-<  ROCKET  PERFORMANCE 

CARD  NO.  2  ;mT0MS  =  AL, H  ,N.CL, 0, 0  SFL0R=F  CCKET»EQ 

CARO  NO.  3  iP  =  -*t  ,  MM  ‘FLAGS,  11  =  50  i  PRINT,  SUM -‘A  cy 

CARD  NO.  t  ;CHAMdfc'R,CREM=AL,MEAr-C,WT  PCT=17, 

CARD  NO.  5  CHt.isM^NCLOAjHEATs-TC  .73,WT  PCt=6l, 

CAPO  NO.  6  CHEM=C.ll’H9.?55  =  31.5S--N.l3ct,HEAT=-62.6,WT  PCT=18 

CARD  NO.  7  i E N 0 


This  input  sample  is  a  wind  tunnel  problem  for  the  calculation  of  equi 
librium  conditions  and  then  the  calculation  of  the  nonequilibrium  chemistry 
after  departure  from  equilibrium.  The  flow  type  is  defined  on  card  2.  The 
$M,  $EQ,  and  $KF  data  streams  are  for  use  in  the  nonequilibrium  calculation 
The  nozzle  shape  is  defined  by  the  $AREA  VS  X  data  stream  starting  on  card 
15.  The  $PRINT  data  stream  is  for  output  during  the  nonequilibrium 
calculation;  the  equilibrium  portion  will  automatically  print  answers  at  48 
stations. 


CARO  no. 
CARO  NO. 
CAPO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARD  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 


a 

9 

10 

11 

12 

13 

IV 

IS 

lb 

17 

ia 


*TITL£=WlNOTUNNEL, 
JFLOW=HlNil,£Q,NGN£ 
JP=1200,ATM  ST  =  ’ 
$M0L/GM,02=. 539335 
N=. 25138C J 
ifl  2=0 

♦  Ml  = 

♦  M2  = 


23. 
23  » 


iM  1=0  2 
i£Q=02, 

$£0-02 i  ... 

§2Q=02  »  ♦  M3  =  23, 
|£Q=N2 ,  ♦  M4  =  — 
$£Q=N0,  ♦  M4  =  N, 
$£Q=0,  ♦  N2  =  N,  ♦ 
i £0=0  2,  ♦  N  =  3 ,  ♦ 
$EQ=02,  ♦  N2  =  2ND 
(AREA  VS  X,  .73  51b 
. 11402,1. 5494.  .11 
.56312,5,  .96224.7 
SPRINT, I NT cfiV Ac =20 


iaulLIdPIUM  A  NO  NON  EQUILIBRIUM 
3  $AT0MS=0»N 
350 

6»£-3.N2  =  .2415C3  04E-l  ,N0=. 3940 631  IE-2, 
5£-2,0=.9537491E-2 
IM  3=  N,N2 , NO  |M4  =  * 


iKF=.1112E26, -2.5.-.1 17204E6 
JKF=.1112E26,-2.5 ,-.l  17  204E6 
SKF=.lll2F26,-2.5,-.ll7204E6 
iKF=.7338E26,-2.5 ,-.22-.475E6 
M4  JKF=.2296 E26 ,-2.5, “.148 98  (E6 
f KF=1.E12, .5,-75487. 

$KF=.  10  14  £12, .5,-6198. 

f<F=  .26  £9El C,-.  106874E6 
,  .265  17,1.  .12438,1.5. 

852, 1.651,  .14872,2,  .26068,2.5,  .262 
,  1.67464,10 

, MAX  ANSMERS=200, XSTOP=10  SENO 


»  Ml 
♦  M2 
f  M3 
►  M4 

*■  3.  ♦ 

MO 

MO 


>20  2,  3, 


The  purpose  of  this  equilibrium  rocket  performance  calculation  is  to 
calculate  the  area  at  the  nozzle  throat.  The  temperature,  pressure,  and  cold 
pressure  must  be  input.  In  addition,  the  mass  flow  and  heat  of  formation  of 
each  system  component  must  be  input  by  means  of  the  $CHAMBER  data  stream. 

It  is  often  of  interest  to  run  the  reverse  of  this  case,  that  is,  to 
input  the  throat  area  and  total  mass  flow  of  the  system  so  that  the  program 
can  calculate  the  chamber  temperature.  If  the  total  mass  flow  and  area  are 
specified  by  means  of  flags  18  and  IS.  respectively,  it  is  assumed  that  the 
chamber  temperature  is  to  be  calculated. 


CARC 

CARD 

CARD 

CARO 

CARO 

CARD 

CARD 


NO. 

NO . 
NO. 
NO. 
NO  . 
NO  . 
NO. 


3 

v 

5 

6 
7 


„  T..Tc  -  =r0C<£T  NDZZ.E  THROAT  Af-fA 
.r_JW  =  RCC<-:T,n.  SA7CNS=H£,  S,F,  0 


I  Ai.  C  UL  AT  :  ON 


jp=cw; , -h,  : old=329  ;  '  iP^iNT,  throat 


ChAM3E= . ,  CHtM=HE  ,  -UiC  F  L  0W=  2  .  2  »  H  £  AT  =  J 


CH£;-1=5F6  i  MASS  r  L  OR  =  2 . 7  ,  H  £A  T=  -Z  9 1 . 7, 


C"tEM=02 . MASS 
id  NO 


F-0W  =  .57  ,HEAT= 


J 


The  Input  for  the  General  Solver  case  and  the  Mollier  case  presented  here 
are  very  similar.  The  General  Solver  case  uses  the  $INCL  data  stream  to  limit 
the  total  number  of  species  in  the  system  to  the  five  species  listed  in  this 
data  stream.  If  the  $INCL  is  not  present,  all  species  on  thermodynamic  tape 
comprised  of  N  and  0  atoms  would  be  included  in  the  system. 

In  the  Mollier  solution  case,  a  fuel  and  an  oxidizer  are  present  with  a 
mixture  ratio  of  28.7.  Note  that  the  label  FUEL  is  included  after  the  methane 
is  defined  so  that  the  program  can  identify  which  of  the  two  species  is  the 
fuel. 


CARO  NO. 
CARO  NO. 
CARO  NO. 
CARO  NO. 
CARD  NO. 
CARD  NO. 
CARO  NO. 


iTITLE  =  GcN£RAL  C0_\/£5 
iATGMi=N)0  i  r  L  Jrf  =  jEN  S  0  L  V  c.  r. 
iT  RANGh=63CC, 25v  ;P  RANGE=12uC»  12Q  C »  5 
*INGL  =  N,C,N2,G2,‘J0 

SMJL/GM, 02=.5  29J:sa7E-3f N2=. 2R 1 5C 3 0% E -1 , 

N  0  =  •  3  9**  I  6  3 1  102,  N  =  .  25193CC5  02,  0  = .  95  3  7C  c1E- 2 
*cNO 


CARC  NO. 
CARO  NO. 
CARD  NO. 
CARO  NO. 
CARO  NO. 


iTITLE=MQi.lIO  SA.-?_£  CASE  •  A  T  C?1S=  C  ,H  ,  N ,  0 
5FL0W  =  M0LLIER  IT  MNGE=2w2C,1S0C,2; 

*PRANGE  =  U  ,9, 1 

lEP.CHEMx";  ~t,,  _  A  3  E  L=  ME  T  H  ANE  ,  R  UEL  ,  MI  X=28 . 7  , 
CHE  M=0«A2u5N  1.58,  -  A 1 E  L=A  If  i£NO 


APPENDIX 


NEST  EXPANDED  VERSION  -  AEROSPACE  ONLY 


This  version  of  NEST  was  developed  for  the  running  of  nonequilibrium 
cases  involving  large  number  reactions  and  species. 

RESTRICTIONS 


1.  2000  total  number  of  reactions 

2.  125  species 

3.  15  elements 

4.  20  catalyst  sets  with  <  100  species  per  set 

5.  One  global  catalyst 

6.  16  flash  lamp  functions 

7 .  A  species  may  appear  in  at  most  250  chemical  reactions 

8.  8  analytic  flash  lamp  functions  are  available 
ADDITIONAL  SPECIAL  OPTION  FLAGS 

16  *  1 ,  the  long  printout  of  the  forward  and  reverse  rate  constants  plus 

the  reaction  species  that  appears  before  the  actual  NEST  calculation 
begins  will  be  suppressed. 

33  x  coordinate  where  rate  printout  is  to  start.  (Omit  RATES  on 
$PRINT  input.) 

34  x  coordinate  where  rate  printout  is  to  cease. 

53  *  1,  suppress  the  printout  for  the  rate  of  production  for  each 

species  at  each  print  station. 

56  *  0,  vibrational  species  (default) 

-  1,  HF(v, J)  v-0,6  J-0,40 

■  2,  DF(v,J)  v-0,6  J«0,40 

*  3,  H2(v,J)  v«0,l  J*»0,15 

■  4,  D2(v,J)  v*0,l  J«0,15 

-  5,  HD(v, J)  v-0,1  J-0,15 

For  a  calculation  using  more  than  one  of  the  above  species,  i.e., 
HF(v, J)  +  H2(v,J),  $FLAG56-13 


101  -  Analytic  flash  lamp  number  (1-8) 

102  *  Input  information  for  analytic  flash  lamp 

150  - 
TIMING 

Note  that  increased  running  time  is  most  dependent  on  the  total  number  of 
species  in  the  system  and  less  dependent  on  the  total  number  of  reactions. 

CONTROL  CARDS 

The  following  files  must  be  attached  if  the  expanded  version  of  NEST  is 

used 

4MAINEXPAND,  ID-06629 
4SUBSEXPAND,  ID-06629 


I 


LABORATORY  OPERATIONS 


I 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting  exper¬ 
imental  and  theoretical  investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyalcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  hlgh-temperature  thermomechanics,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chalcal  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  mil lime ter- wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  Forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  Mterlals  in 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 
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